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In ref. [T] we have used the Zagreb realization of Carnegie-Melon-Berkeley coupled-channel, uni- 
tary model as a tool for extracting pole positions from the world collection of partial wave data, 
with the aim of eliminating model dependence in pole-search procedures. In order that the method 
is sensible, we in this paper discuss the stability of the method with respect to the strong variation 
of different model ingredients. We show that the Zagreb CMB procedure is very stable with strong 
variation of the model assumptions, and that it can reliably predict the pole positions of the fitted 
partial wave amplitudes. 

PACS numbers: 14.20. Gk, 12.38.-t, 13.75.-n, 25.80.Ek, 13.85.Fb, 14.40.Aq 

I. INTRODUCTION 

Various analyses of meson-baryon scattering designed with the goal to extract resonance properties from 
data are available in literature. They fundamentally differ in technical approach how the model is imple- 
mented (effective Lagrangian approach, K matrix formalism, K matrix approximation, phenomenological 
coupled-channel formalism, t-channel dispersion relations, hyperbolic dispersion relations, e.t.c), however 
most of them strongly insist in obeying fundamental physics principles like Lorentz invariance, crossing 
symmetry, unitarity and analyticity. In ref. PQ, we have extracted the pole positions from a "world collection" 
of partial wave data and partial-wave amplitudes understanding them as partial-wave data, with the aim 
of eliminating the model dependence from the individual pole-search procedures. We have analyzed some 
major PWA amplitudes like: i) Karlsruhe-Helsinki (KH80) amplitudes [2! for the 7rN elastic scattering, 
where fixed t-channel dispersion relations are combined with Pietarinen expansion to impose and maintain 
point-to-point analyticity; ii) VPI/GWU single energy (GWU-SES) [3], and energy dependent solutions (WI08) 
[3j |4] where Chew-Mandelstam K-matrix formalism together with dispersion relation constraints is applied; 
iii) dynamic coupled-channel Lagrangian approaches from EBAC model [SJ [5] which is based on an energy 
independent Hamiltonian which is derived from a set of Lagrangians by using a unitary transformation method 
[7j; iv) Julich model [Hj which is also a dynamical coupled-channels model of meson production reactions 
in the nucleon resonance region which includes 7rN, 77N, 7rA, ctN, and pN channels; v) Dubna-Mainz-Taipei 
(DMT) model [91 UOj . a meson-exchange model for pion- nucleon scattering which was constructed using a 
three-dimensional reduction scheme of the Bethe-Salpeter equation for a model Lagrangian involving it, 77, N, 
A, p, and a fields; and vi) and the Giessen K-matrix approximation (analyticity violating) model |11) . Let 
us not fail to mention other, important theoretical approaches like Cutkosky CMU-LBL phenomenological 
coupled channel model [T5], Pittsburgh phenomenological coupled-channel model [T3] and Bonn-Gatchina 
K-matrix approach without dispersive part |14) , which are for different reasons not included in analysis of ref. [T] . 

In this paper we discuss how stable our recommended procedure is with respect to our model assumptions 
in order to give an quantitative estimate how confident the results of that report actually are. 

We have been triggered to start the whole enterprise by the present uncertainty how to make the least model 
dependent comparison between scattering theories which analyze experimental data on one side, and QCD 
models or lattice QCD calculations on the other. In other words, we have been dealing with the problem how 
to precisely define what a resonance actually is. And it has to be done on both sides, on scattering theory side, 
as well as on the QCD side. 

One, and the most intuitive way to define a resonance in a scattering process, is to define it as an 
intermediate state of two particles when they, because of attractive interaction, dwell in the vicinity of each 
other longer then in a standard scattering process. When this definition is transformed into the language 
of scattering theory, since 1976 it seems that, as Holder has said in [IS]: "It is 'noncontroversial among 
theorists' (see Chew 176) / and the references in my 'pole-emics', p. 697 in ref. \17j ) that in S-matrix theory the 
effects of resonances follow from first order poles in the 2nd sheet. " However, if resonances are to be identified 
with the first order poles in the 2nd sheet, we can not but observe that characterizing such states, namely 
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identifying and quantifying scattering matrix poles in the complex energy plane turns out to be quite a challenge. 

Different methods aimed onto identifying one of the possible repercussions of a resonant state onto 
measurable quantities have been developed. First it has been speculated that a Brcit-Wigner function 
should be a good representation of the complex energy pole, so the Breit-Wigner formula for spin zero 
particles together with its generalization to the non-zero spin case has been developed in 1936 (see an 
illustrative discussion in Cottingham and Greenwood, p. 241 in ref. [IB])- Then in [T^ Dalitz and Moorehouse 
considered the scattering matrix eigenphases, and their rapid increase through tt/2 have been taken as 
a signature of pinpointing a scattering matrix pole (namely, the corresponding eigenvalue of the reaction 
matrix K = i (S-1)/(S+1) then has a pole at this energy [2Q]). Finally, closely related to the afore two, the 
quick backwards looping of Argand diagram have been also extensively discussed by Dalitz and Moorehouse [T5] . 

As the time passed by, the real meaning of these methods, namely the fact that they are devised only 
to reach the scattering matrix pole in the complex energy plane in a certain approximation, tended to be 
overlooked. These methods alone have been identified as resonance definition procedures, and what they found 
has been proclaimed to be a resonant state. However, their inherent model dependence, like for instance 
unitary addition of non-resonant background, has been systematically neglected (see the elaboration on 
Breit-Wigner model dependence in ref. [17 ). The problems of eigenphases, linked to the manifest violation of 
Neumann- Wigner no-crossing theorem |21) when rapidly going through 7r/2 even for the case of simply constant 
background terms, have been suppressed in spite of the extensive elaboration of the problem given by Dalitz and 
Moorhouse in ref. [19j. They have explicitly said that: "With such a complexity of branch cuts without physical 
significance, we must conclude that the eigenphase representation for the S matrix is not generally a useful 
representation for the scattering in the neighbourhood of a resonance", but have anyhow allowed the possibility 
to use it in a more restricted situation: "However, we should emphasize also that the eigenphase representation 
is a perfectly acceptable and economical representation of the S matrix in situations where the eigenphases 
vary slowly and do not cross". The single-channel nature of pole search methods (with the exception of 
the time delay which is in principle of multi-channel character) also did not represent much of a prob- 
lem, because experimental data in inelastic channels were scarce. The fact that a single- channel procedure can 
not reveal much information about resonant states being 'far away' from that channel, has also been suppressed. 

In ref. pQ we have returned to the root of the problem - we have fully respected that we are looking for 
the pole in the complex energy plane having at our disposal only data on the real axes, and instead of using 
different pole-extraction methods to extract a pole per resonance, we have suggested to use only one, Zagreb 
CMB method [251 US] an d extract poles from a selected collection of partial wave amplitudes simultaneously. 
We have decided to use only one method to extract pole positions from all published partial waves analyses, 
and inspect the result with the aim to distinguish which part of the disagreement in pole positions is coming 
from the different analytic structure, and which is coming from the different input (different PWA functions). 
In other words, we have taken all sets of partial-wave amplitudes, treated them as nothing else but a good, 
energy dependent representations of all analyzed experimental data, and extracted the poles which are needed 
by the CMB method. We observe that even when practically all analyzed PWA use very similar input data 
set, their PWA solutions do differ in spite of reporting the similar quality of fit to the input data (similar 
reduced x 2 ). So, from the theoretical point of view, they equally well describe the experiment. And now 
these, however similar, but still different curves through different analytic continuations of each model generate 
corresponding and different sets of poles. So, in addition to the issue of slightly different input, the error of 
unknown analyticity is superimposed to it. It is important to realize that these poles, even for the identical 
set of input amplitudes, should not necessarily coincide, because the models used for analytic continuation 
are intrinsically different in their analytic form (for illustration see ref. |24j ) . The idea of using only one 
model (Zagreb CMB) to extract the set of poles from different PWA treating them as partial wave data boils 
down to testing the internal agreement of input data sets. In this way, the difference between poles of various 
solutions is attributed only to the under-determinacy of input data and not to the analytic structure of the 
models in question. Simply, different poles obtained in this way quantify the difference in PWA solutions with 
respect to the similar input, and disregard the different analytic form used to obtain them. In this manner, 
all uncertainties originating from different analytic properties of different models are avoided, and the only 
remaining errors are the quality of the input and the precision of CMB method itself. Therefore, averaging 
and error analysis of pole positions is sensible, and can be safely carried out. To answer the question of a 
correct choice of analytic form is a more complex problem and will be addressed elsewhere. 

In this paper we discuss how sensitive this method is to CMB model assumptions when our input data set is 
well defined. We have first chosen the reliable input data set, a data set for which the minimization is fast, stable 
and unambiguous. Then we have identified our model assumptions, specific ones which we may check, and more 
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general ones which are unfortunately beyond our reach. The general ones, like the isobar character of the model 
representing the unstable channels with a quasi two body channel with decay properly taken into account, could 
not be checked directly. Other ones, like for instance the form of the channel propagator (channel-resonance 
vertex function), dependence on the dispersion integral subtraction constant, background parameterization, 
number of channels, mass of the effective channel, etc. we check directly. Instead of taking different recipes 
for each model assumption, we have simply modified the used functional form with an completely arbitrary 
function which is just changing the range of interest, and repeated the fit. In this way the stability with respect 
to different model assumption is performed. 

II. GENERAL IDEA 

If we want to apply a CMB formalism with the intention of extracting pole parameters from a "world 
collection" of partial wave data and partial wave amplitudes, one has to answer a natural question how stable 
the formalism itself is with respect to model assumptions. We are aware that we can test the sensitivity of 
the model upon its ingredients, but the reliability of the isobar model itself can not be tested in this way. 
Therefore, at this instant, we still can not say how strongly our answer depends on the fact that we assume 
that the interaction mechanism of isobar intermediate states is the only relevant one. 

Our present analysis involves yet another approximation, and that is that possible three body states are 
effectively represented by a group of two body ones. We shall not test this hypothesis either, even knowing 
that according to some considerations these effects can build up to 30%. 

Even before starting the procedure, we have to answer a nontrivial question: 
"What is the optimal input data base so that our pole positions are unique for the model without any 
modifications? " 

Once this answer is given, we have to analyze how the final pole positions depend on model assumptions like 
(i) shape of the channel-resonance form factor; (ii) choice of the dispersion integral subtraction constant; (iii) 
masses of the effective channel; (iv) number of channels; (v) type of the background parameterization. 
The driving idea of this paper is: 

"To drastically modify only one part of the model at a time, and then make a fit to the chosen data base. The 
dissipation of pole positions known from the original fit will give us information how important certain part of 
the model is for the stability of the solution. " 

We will just add that these changes are intentionally artificial, aimed only at obtaining the effect, and not 
paying attention to the physical meaning. 

III. FORMALISM 

The CMB model is isobar, coupled-channel, analytic, and unitary model, where the T matrix in a given 
channel is assumed to be a sum over the contributions from a number of intermediate particles (resonance 
and background contributions). The coupling of the channel asymptotic states to these intermediate particles 
determines the imaginary part of the channel function, and is represented effectively with a separable function. 
The real part of the channel function is calculated by the dispersion relation technique, thus ensuring analyticity. 
Besides the known resonance contributions, the background contributions are included via additional terms with 
poles below the 7rN threshold. Due to the clear analytic and separable structure of the model, finding the pole 
positions in CMB model is trimmed down to the generalization of the dispersion integral for the channel 
propagator from real axes to the full complex energy plane, and this is a very well defined procedure. In 
practice, we instead use a very stable, and numerically much faster analytic continuation method based on the 
Pietarienen expansion |25] in order to extrapolate the real valued channel propagator into the complex energy 
plane. 

A. Formulae 

Our current partial-wave analysis |22j is based on the manifestly unitary, multichannel CMB approach 
of ref. |12j . The most prominent property of this approach is analyticity of partial waves with respect to 
Mandelstam s variable. In every discussion of partial- wave poles, analyticity plays a crucial role since the poles 
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are situated in a complex plane, away from physical region, and our measuring abilities are restricted to the 
real energy axis only. To gain any knowledge about the nature of partial-wave singularities would be impossible 
if partial waves were not analytic. Therefore, the ability to calculate pole positions is not just a benefit of the 
CMB model's analyticity but also a necessity for resonance extraction. In this approach, the resonance itself is 
considered to exist if there is an associated partial- wave pole in the "unphysical" sheet. 



We use the multichannel T matrix related to the scattering matrix S as: 

S a b{s) = 5 ab + 2iT ab (s) 1 

where 5 ab is Kronecker delta symbol. The T-matrix matrix element is in the CMB model given as: 

N .IL 

T J ab L { S ) = £ fl L {s)^pJs)lttG J l3 L { S )l J 0b L yfp&)ti L {8) (1) 

where a(b) represents the outgoing (incoming) channel. In our analyses we use a,b = irN, r]N, ir 2 N. The 
initial and final channel b(a) couple through intermediate particles labeled i and j. The factors ji a are energy- 
independent parameters occurring graphically at the vertex between channel a and intermediate particle i and 
are determined by fitting procedure. Also occurring at each initial or final vertex is form factor f^ L (s) 

fa L (s) = (— %=~ ) , (2) 



Qla+VQla+Ql, 



and phase-space factor p a (s) 



Pa{s) = — 7r , (3) 

V s 

where s = W 2 is a Mandelstam variable, and g a (s) is the meson momentum for any of the three channels given 
as 

_ y/js - (m + m a ) 2 )js - (to - m a ) 2 ) 

Qa{S)- • W 

Furthermore, L is the angular momentum in channel a, and Qi a , Q2a are constants. The factor /^(s) 
provides appropriate threshold behavior on the right-hand cut, and also produces a left-hand branch cut in 
the s plane. Parameters Q\ a and Qia are chosen to determine the branch point and strength of the left-hand 
branch cut. In our analyses they have been taken to be the same, and are fixed to the mass of the channel 
meson a. 



G^j is a dressed propagator for partial wave JL and particles i and j, and may be written in terms of a 
diagonal bare propagator G^ L and a self-energy matrix S^ ; L using Dyson equation 

N .JL 

G{ 3 L is) = G^ L is)+J2 G$Z L (s)Ztf(s)Gtf(s). (5) 

k,l=l 



The bare propagator 

g ojl {s) = (6) 

Si s 

has a pole at the real value s*. The sign a = ±1 must be chosen to be positive for poles above the elastic 
threshold which correspond to resonance. 

The nonresonant background is described by a meromorphic function, in most of the cases consisting of two 
terms of the form ^ with pole positions below ttN threshold. For that case, the signs of the terms are opposite. 
The positive sign correspond to the repulsive and the negative sign to the attractive potential. In principle the 
number of poles can be increased arbitrarily (see the next subsection on background representation), but in 
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reality the number is never larger than three. 

£^ is the self-energy term for the particle propagator 

z{t(s) = J2rla-<s>i L (s)-jL L 



(7) 



The $^ L (s) are called "channel propagators". They are constructed in an approximation that treats each 
channel as containing only two particles. We require that have, in all channels, correct unitarity and 
analyticity properties consistent with a quasi-two-body approximation. 



The imaginary part of (s) is the effective phase-space factor for channel a: 

Irakis) = [f a L {s)fp a {s). 
The channel propagator is evaluated on the real axes only 

i2L+l 



Im<I>(:c) 



[q(x)Y 



2L ■ 



(8) 



(9) 



where by x we stress the fact that values are on the real axes. The real part of (x) is calculated using a 
subtracted dispersion relation 



Re$(x) = Re$(x ) 



x — Xq 



\m®{x')dx' 

{x 1 — x)(x' — Xq) 



(10) 



where x a = (m + m a ) 2 and Re$(xo) is a subtraction constant at the subtraction point xq. In all publications 
[2H [53] we have chosen the convention that x = x a and Re$(x ) = 0, but this is not necessarily so. For 
better understanding, the structure of the channel-intermediate particle form factor and its relation to the full 
T-matrix is given in Fig. [T] 

T ab a 
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FIG. 1: (Color online) Parameterization of channel-intermediate particle vertex in CMB model. 
We give a matrix form of the final T-matrix as defined in Eq. ([I]) : 

G (s) 



T(s) = A/im$(s) • f 



I- 7-$(s)-7T .g o(s ) 



7 • y Im<f>(s) 



(11) 



B. Background contributions 



Theoretically, background term is in principle a free function of energy. However, an approximation is very 
often used that one may decompose the background term into a finite number of unphysical pole terms, i.e., the 
background is a meromorphic function. This approach is very convenient for the CMB type calculation because 
integral equations for the background poles are solved in the identical manner as for physical ones, so has it 
been applied and strongly defended by Cutkosky et al [T2] in CMB PWA. Zagreb group has taken it over in [22] , 
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Let us elaborate on what we have actually done when we assumed that the background indeed is a 
meromorphic function. In that case, instead of solving a nonlinear equation for the T-matrix pole positions 
when the background has a general energy dependent form, we have actually introduced a new type of poles 
which mimic the background. Here we rely on the approximation that the solutions of a nonlinear equation 
with general type background will be close to the solutions of "pole-type" equation, the equation when the 
meromorphic background is very close to the original function. Consequently, we do not have one, but two 
types of T-matrix poles: i) poles which correspond to the physical internal singularities; and ii) poles which 
correspond to the background. We call the first type genuine poles, and the second ones we call dynamic poles. 

The speculative question is: "How realistic the meromorphic background approximation actually is?" 



C. Pole extraction procedure 

The poles of the T matrix given in Eq. are found by solving the equation: 

det G _1 (s) = det [e y %(sj - s) - S y ] = 0. (12) 

This ia a complex equation, and can be solved only by analytic continuation of the physical, experimentally 
accessible T-matrix values which of course lie on the real energy axes into the complex energy plane. Let us 



observe that the only part of the Eq. ( 12 ) which is not defined in the complex energy plane is the channel 



propagator which is defined on the real axes only. 

Therefore, the problem of finding the poles of Eq. ([T]) is reduced to analytic continuation of the channel 
propagator. 

This analytic continuation can be done either by numeric principle value integration, or by Pietarinen ex- 
pansion. The numeric integration of the channel propagator, used in ref. [22] , is a straightforward procedure, 
but requires a lot of computer time. Instead, we have used a more elegant way, a Pietarinen expansion used by 
Karlsruhe-Helsinki group [2J. 

From Eq. (9| it is evident that $(s) has a square-root type singularity in cm momentum, and analytic contin- 
uation into the complex energy plane should take it into account. Instead of calculating the dispersion integral 
( 10 1 for each point in complex plane, we decided to expend the function $(s) in power series of a new variable, 
which accounts for all analyticity requirements. We use the expansion (similar to Pietarinen's in Ref. |25| or 
Ciuli's HH) 



$ I ( s ) = J2Cn [zm\\ (is) 



n=0 

where C n are coefficients of expansion. The new channel dependent variable is given by its principal branch 



Zl {s) = VIM- ( 14 ) 

a + y/X a - s 

with the tuning parameter a. This function is fitted to a data set consisting of imaginary parts of $(x) from 
Eq. ([9| and real parts of $(x) calculated from dispersion relation (10 1, both of them evaluated at real axis 
(hence x). 

The general idea is that the $(s) inherits analytic structure from Z(s). We obtained parameters a and 
coefficients C n for each channel, and for all analyzed partial waves. The least-square fit is considered to be 
good if it meets following conditions: (i) small number of coefficients C n needed (7 or 8, at most), (ii) the 
function fitted to the part of data set, when extrapolated outside of the fitted region is consistent with the 
rest of data, and (iii) fitting just imaginary part of <i>(a;) produces real part that is in agreement with values 



obtained from ( 10 1 



The channel propagator given by expansion (131 is obtained quite accurately on the real axes, and as the 



expansion possess a correct analytic structure of 2— >2 scattering processes, the deviation between expansion 



(131 and the correct value given by Eq. (10) should be fairly small for the resonant region in the vicinity of 
physical axis. 
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Every channel opening is responsible for two Riemann sheets: the first (physical) sheet with physical partial 
waves, and the secondary (unphysical) sheet with resonant poles. To get to the unphysical sheet it is enough 
to use the second branch of Z(z) 



Z n (s) = (15) 

a - v^q ~ s 

Finally, it is evident from Eq. ([I]) that all poles of each partial wave must by construction be the same in all 
channels and, in fact, equal to the poles of the resolvent G(s). 

In practice, for the S-wave 7rN channel propagator, the realization of the Pietarinen expansion on the real 
axes is shown in Fig. [2| Red points represent the imaginary part of the channel propagator $ on the real axes, 
and are given by Eq. (|9|), blue points represent the real part of the channel propagator $ on the real axes, and 



are calculated by the dispersion relation (10), and the full line represent the function $i(x), a complex function 
defined on the full complex energy plane, and having a proper branch points and corresponding cuts for each 
channel. 




FIG. 2: (Color online) Pietarinen expansion for the Sn partial wave. 



As the channel propagator is now known in the whole complex energy plane, the equation Eq. (12 1 was 
solved num erically as the complex function of the complex arguments. In reality, instead of looking for zeroes 
of Eq. (12 1 we have looked for the poles of det G(s) _1 . In practice we have calculated |detG| as a function of 
Re(s) and Im(s) - (3D plot) specifying exactly the Riemann sheet in question (we have defined the appropriate 
momentum square roots in the channel propagator). Then we have looked for the values of Real(s) and Imag(s) 
when |detG| becomes bigger then an arbitrary number. 

We regarded the obtained values of Re(so) and Im(so) as pole positions, and illustrate the procedure in Fig. [3] 

\det G\ 




FIG. 3: (Color online) The principle of numerical finding the pole positions using Eq. ( 12 1 with the Pietarinen expansion 
of the channel propagator. 
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IV. SOURCES OF MODEL DEPENDENCE 

Understanding that the general assumptions of the model like its isobar, two body character in which the 
three and more body channels are identified with a quasi two body channels are properly taken into account, we 
now identify other sources of model dependence which at the present moment represent additional, systematic 
and uncontrollable uncertainties: 

1. the form of the channel propagator (meson- resonance vertex function) 

• its asymptotic part 

• its intermediate part 

• threshold behavior 

2. dispersion integral subtraction constant 

3. background parameterization 

4. number of channels 

5. mass of the effective channel 

The idea of present work is to allow for the extensive change of the chosen model features, and see at which 
level the pole positions are starting to move. The extensive change of the chosen model feature is realized by 
multiplying it with a freely chosen function (no physical meaning whatsoever) which influences only the tested 
part of the model assumption, and parameters are varied in such a way that the change in the analyzed part 
in/decreases from ~ 10% to more than 100% of the original value. 

Prerequisite of such a procedure is that the input data base is such that it ensures a confident and reliable 
fit. Namely, we have to be positive that the variation of the extracting pole positions is coming from changing 
the ingredients of the model drastically, and that the fit to the input data base converges fast and uniquely. 
Therefore, we first describe the choice of the input data base, and then we test the ingredients of the model. 

V. TESTING THE MODEL DEPENDENCE 

A. Defining the input data set 

If we want to test the stability of the model with respect to its ingredients, we first have to choose such an 
input data base (set of input data our modified model will fit) for which the fitting procedure with original 
model is fast, unique and reproducible. In other words, we want to choose such an input data base for which 
the original model reproduces its own, input values safely, quickly and unambiguously. In that case, we can 
vary the ingredients of the model, and claim that the resulting shift in pole positions is the result of change 
in the model itself, and not the consequence of inadequate input data. This is in a way a test of internal 
consistency of the model, the proof that the model when fitting its own result obtains the parameters that 
were used to generate the input (logistic tautology). However, when we start changing the ingredients of the 
model, the analytic structure of the model itself is changed, so we necessarily expect to obtain a different set 
of resulting pole parameters. The deviations of these parameters from the original input values are now the 
measure of sensitivity of the model to different model assumptions. 

As the input for this work we have chosen partial wave amplitudes of Zagreb model [22] . Once the input 
parameters are chosen and fixed (bare poles and channel-resonance coupling constants 7), we have all partial 
waves at our disposal. In this section we test how stable the fit is for the first three lowest angular momentum 
partial waves Sn, Pn and D13 with respect to how many channel data sets we have used. 

In Zagreb model we have three set of partial wave amplitudes at our disposal: ttN — > ttN, ttN — » tjN and 
ttN —> ir 2 N. In principle, all Zagreb PWA are smooth energy dependent functions. We could have fitted 
these solutions directly. However, we have been of an opinion that by doing so it would introduce too many 
analytic constraints coming from the fact that these solutions "memorize" the analytic form of CMB formalism 
used to obtain them. We wanted to eliminate these additional constraints, so we have distributed the energy 
continuous Zagreb solutions by the recipe given in ref. [27]. 
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Instead of using smooth theoretical curves and avoid pathologically small Xr > we constructed pseudo data 
points by normally distributing the model input in order to simulate the statistical nature of really measured 
data. The standard deviation a was set to 0.02, similar to the average error value of GWU data [31 H|. 
By using this procedure we were able to produce a set of ttN — > irN, irN — > r/N and ttN — > tt 2 N partial-wave 
data which reproduced the experiment, and when fitted, they gave realistic \ 2 values comparable with those 
obtained by GWU/VPI SES fits. The pole parameters of the input data set are given in Table [T] The input 
data sets together with our fit are shown in Fig. |4| 

To simplify the discussion, we introduce the generic abbreviation L2/2j(mass) for each pole position produced 
in Zagreb model. 

To illustrate, the Sn (1518) pole denotes the first, 1518 MeV Sn pole, and is identified with the N(1535) 
Sn PDG resonance. Observe that all Zagreb solutions correspond to an existing PDG resonant state, with the 
exception of two Pn intermediate energy solutions. In Zagreb model we get two Pn solutions: Pu (1708) pole 
which is in [2"2l |2"5] identified with N(1710) Pu state, and a new solution P n (1728) which is not confirmed by 
PDG. 

In Tables [TTJ |III[ |IV| we show the results of the test in which we have looked for the optimal number 
of fitted channels required for the stable fit. Results obtained by fitting only elastic channel are denoted 
by Fit 1. Fit 2 stands for a case when the ttN — > r/N channel is included in a fit in addition to clastic 
channel. Results denoted by Standard fit are obtained by fitting all three channel simultaneously. The 
reduced chi-squared for the limited data input (one channel for Fit 1 and two channels for Fit 2) is denoted 
with \ 2 - We have calculated the new chi-square parameter which is obtained when solutions for Fit 1 and 
Fit 2 are used to calculate the agreement with all three channels without any fitting of the third channel, 
and denoted it as tot x\- This new parameter now describes how well single and two channel fits reproduce 
all channels at the same time. We are aware that fitting two out of three channels should in principle 
fully predict the third one for the manifestly unitary model as Zagreb CMB, but this is true for the ideal, 
theoretical and smooth input. However, when we have decided to distribute the smooth input by adding 
realistic errors, in spite of the fact that the statement is still in principle valid, the reality makes the fit rather 
complicated and slightly unstable if only two channels are fitted. That can be seen when analyzing the tables 
containing the \ 2 an d tot x\ f° r the two and three channel fits. Both parameters are qualitatively similar, 
but for the three channel fits, \ 2 is still always slightly better then the tot X 2 R parameter for the two channel ones. 

Hence, we have shown that reduced \ 2 indeed is sufficiently stable if we fit two out of three channels, but 
in spite of the fact that fitting three out of three channels effectively is over-constraining the problem, we have 
decided to use it thought this paper because this technically simplifies the fitting procedure when generated 
pseudo-data are fitted, and also opens the possibility to increase the number of channels. 
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FIG. 4: (Color online) The input data sets together with our fit to them with the standard model are shown for the Sii, Pn and 
D13 partial waves T-matrices. 



TABLE I: The original parameters of the Zagreb solution for the Sn, Pn and D13 partial waves. 



Partial 


Bare poles 


Dressed poles 


waves 


W Sl W S2 W S3 W S4 
MeV 


( RcW \ 
l v -2ImW J 


( RoW \ ( 
l v -2ImW y ) l v - 

MeV 


RoW \ 
-2ImW ) 


( RcW \ 
l v -2ImWy 


Sn 


1523 1640 1837 


(■so 




^1784\ 
v 424,1 




P11 


1607 1771 2186 2852 


f 1360\ 
{l6l) 




f 1728\ 
v 138 j 


/ 2114\ 
V 345 ) 


D13 


1582 1880 2497 


(■50 




f 1942\ 
V 47lj 


/2699\ 
^ 558 ) 



TABLE II: The stability of the fitting result for the Sn partial wave. 





O 


*0 a 




Bare poles 


Dressed poles 








Number 
channels 


Number 
resonanc 


Fitted 
channels 


W S1 W S2 W S3 
MeV 


/ RoW 
l v -2ImW 


\ ( RoW \ / 
j \-2l-mW) \ 

MeV 


RoW \ 
-2ImW J 


Xr 


tot, 2 

Xr 


Fit 1 


3 


3 


1 


1455 1654 2860 


/ 1562\ 

100 ; 


/ 1684\ 
\179) 


f 4090\ 
^2450 J 


1.094 


117.186 


Fit 2 


3 


3 


2 


1522 1635 1823 


{185 ) 


/1640\ 
^ 207 J 


/1790\ 
I 435; 


1.004 


1.711 


Standard 
{it 


3 


3 


3 


1523 1640 1834 




/1642\ 
^ 208 J 


/1779\ 
^ 427 J 


1.016 


1.016 



TABLE III: The stability of the fitting result for the Pn partial wave. 








«" to 

O CD 




Bare poles 


Dressed poles 








Number 
channels 


Number 
resonanc 


Fitted 
channels 


W S1 W S2 W S3 W S4 
MeV 


f RoW \ 
\-2lmWj 


/ RcW \ / 
l v -2ImW y i ^- 

MeV 


RoW \ 
-2ImW^ 


/ RoW \ 
l v -2ImW y i 


Xr 


tot, 2 

Xr 


Fit 1 


3 


4 


1 


1712 1781 2172 3162 


A416\ 
I 69 J 


(T) 


fl788\ 
yl28) 


/2118\ 
\32l) 


0.872 


99.866 


Fit 2 


3 


4 


2 


1734 1923 2191 2330 


f 1360\ 

V 169; 


/ 1695\ 
V 216 / 


fl745\ 
{112 ) 


/2120\ 
^ 350 J 


1.048 


1.739 


Standard 
fit 




4 




1607 1772 2182 2841 








(»«) 


0.958 


0.958 



TABLE IV: The stability of the fitting result for the D13 partial wave. 





O 


O CD 




Bare poles 


Dressed poles 








Number 
channels 


Number 
resonanc 


Fitted 
channels 


W S1 W S2 W S3 
MeV 


/ RoW \ / RoW \ / RoW \ / RoW \ 
l v -2ImW; l v -2ImW y i l v -2ImW y i l v -2ImW y i 

MeV 


Xr 


tot, 2 

Xr 


Fit 1 


3 


3 


1 


1589 1998 2537 


/l500\ /1765\ /1980\ /2720\ 
^ 127 J \ 200 J \ 232 J \ 597 J 


1.295 


137.071 


Fit 2 


3 


3 


2 


1589 1880 2557 


/1507\ A800\ A932\ /2700\ 
^ 127 J \ 126 J \ 508 J \ 595 J 


1.070 


1.555 


Standard 
fit 


3 


3 


3 


1582 1880 2499 


/1506\ /1807\ /1939\ /269l\ 
^ 121 J V 127 / V 485 / 1,583,) 


1.027 


1.027 
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B. The form of the channel propagator $(s) 



The major assumption of the CMB model is done about the functional form of the channel propagator. We 
first define the channel-resonance vertex function in Eq. and then combine it with the phase space factors 
into the analytic channel propagator, see Eqj9|. We take over a recipe extensively used by CMB group [T^] . 
and the full definitions are given by Eqs. |2j |9[and [lO} . 

However, it is very difficult to justify the correctness of these assumptions. So, instead of defending its 
physical reasonability, and discussing the dependence of the vertex function upon the channel resonance 
vertex constants Q\ a and Qi a -, we shall test how sensitive our model is substantial but artificial changes of its 
constituent parts: asymptotic part, intermediate part, threshold behavior. 
Procedure 

(i) We have chosen the modification function completely artificially, requiring that the change of the channel 
propagator is localized only to its constituent part, and that the transition from the modified to non-modified 
part is smooth. 

(ii) The imaginary part of the channel propagator is evaluated. 

(in) The real part of the channel propagator is calculated using the dispersion relation given by Eq. (10 1. 
(iv) The input data set is refitted with modified channel propagator. 



Important: 

We have modified all three channel propagators simultaneously using the same modification junction. 



1. Asymptotic part 



The imaginary part of the channel propagator was kept unchanged at s below 10 GeV 2 . Above 10 GeV 2 the 
correction function g(s) which changes the asymptotic behavior of the channel propagator is introduced: 

/(sH/o + /i-e- Ms - So > 2 , (16) 
Modified imaginary part of the channel propagator (omitting channel indices) now reads: 

T .T, _ flm*(s) for s< 10 GoV 2 /, ? x 

imv„ r - | Imi ( s ).j( s ) for s >io Gov 2 \ L( ) 

Parameters /o, fi and A are determined in such a way to ensure continuity of channel propagator and its first 
derivative at s = 10 GeV 2 , and s — 10.5 GeV 2 . Imaginary and real parts of channel propagator are shown in 
the Fig. §5§ for several values of parameter /o . 




FIG. 5: (Color online) Real and imaginary part of channel propagator. Several values of parameter /o are used: 0.4, 0.7 and 0.9. 



Let us briefly discuss the analyticity and unitarity of afore used piecewise defined functions. When analyticity 
of the channel propagator is imposed via dispersion relations, the analyticity of real and imaginary parts 
themselves is not required [55] • Piecewise defined imaginary part of channel propagator is by construction 
a continuous function having continuous first derivative and not necessarily analytic, but when we use 
dispersion relations to obtain the corresponding real part, we automatically construct the analytic function 
having the same cut. The idea of the paper is to analyze how this change influences pole position of the T-matrix. 
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FIG. 6: (Color online) Real and imaginary parts of Su, Pn and D13 elastic channel propagator. We change intermediate part of 
channel propagator. 

Such a change does not influence unitarity, because in ref. [T^] it has been explicitly shown that the CMB 
formalism by the very construction retains unitary if all used functions are analytic. 

Result: 



As it is seen in Fig. ^ neither imaginary nor real part of the channel propagator are changed in the relevant 
energy range s < 10 GeV 2 . 

Conclusion: 
Pole positions are unaffected. 

2. Intermediate part 

As the intermediate part of the channel propagator, we define the energy domain which is notably above 
threshold, but not exceeding the high energy part above 10 GeV 2 . In this energy range, 1 GeV 2 < s < 10 GeV 2 , 
most of the known nucleon resonances are situated. 



In this part of our analysis channel propagator is multiplied by a function 



9(*) = 



1 i i ?-i 2 

1 + a ■ s + s 

1 + s 2 



(18) 



which changes a shape of channel propagator without changing its asymptotic or threshold behavior. Depending 
on the size and the sign of the free parameter a, the value of the channel propagator in the intermediate range 
2 GeV 2 < s < 5 GeV 2 can be enhanced, or reduced for as much as 100%. 

The imaginary part of channel propagator is given by 

Im $corr = Im $ • g(s), (19) 
and the real part is calculated using the dispersion relation ( |10[ ). 

The change of the channel propagator for the Su, Pn and D13 partial wave is given in Fig. 
Result: 



The shift in pole positions for partial waves Su, Pn and D13 are presented in Tables |V - VII and depicted in 
Fig.Q. 



Conclusions: 



As seen in Fig. ([6| changing the a parameter in ad hoc modification function given by Eq. ( 18 1 within the 
range -1 < a < 1 changes the imaginary part of the channel propagator for up to 100% in the intermediate 
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TABLE V: The extracted Sn partial wave T-matrix poles using different solutions of channel propagator. We change intermediate 
part of channel propagator and the modification is the same for all three used channels 



Solutions 



Bare poles 

W S1 W S2 W S3 
MeV 



Dressed poles 
!!. W ( ReW \ 



-2ImWi 



-2ImWy 

MeV 



ReW 
-2ImW 



Xr 



Sol a = -1.0 
Sol a = -0.5 
Sol a = -0.1 



1522 1640 1837 

1523 1641 1835 
1523 1640 1834 



1505 
176 

1512 
181 

1516 
188 



1636 
208 
1640 
208 
1642 
207 



1773N 
394 ) 
1776\ 
418 ) 
1779\ 
424 J 



1.198 
1.048 
1.031 



Standard fit 



1523 1640 1834 



/1518\ /l642\ ( 

^ 189 y ^ 208 ) \ 



T779\ 
427 J 



1.016 



Sol a = 0.1 
Sol a = 0.5 
Sol a = 1.0 



1523 1640 1834 
1523 1639 1833 
1523 1638 1831 



1518 
189 
1522 : 
190 
/1524 
I 187 



/1522\ /1644\ ( 

\ 190 ) \ 207 ) \ 



1642 
207 
1644 
207 
1645 
210 



1780\ 
427 ) 
1782\ 
433 ) 
1783\ 
445 ) 



1.037 
1.062 
1.105 



TABLE VI: The extracted Pn partial wave T-matrix poles using different solutions of channel propagator. We change intermediate 
part of channel propagator and the modification is the same for all three used channels 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W S3 W S4 
MeV 


/ ReW \ 
\-2lmW) 


/ ReW \ ( 
\-2\mW) \- 

MeV 


ReW \ 
-2ImW j 


( ReW \ 
l v -2ImWy 


Xr 


Sol a = -1.0 
Sol a = -0.5 
Sol a = -0.1 


1603 1772 2181 2929 
1606 1772 2182 2874 
1606 1772 2182 2849 


( 1367\ 

V 139 / 

( 1367\ 
/ 1359\ 


(™) 

(\ 7 7°4 8 ) 


( 1736\ 
V 13 9y ) 
f 1733\ 
V 137,1 
( 173l\ 


/2105\ 
^ 304 ) 
?2111\ 
^ 328 ) 
?2115\ 
V 345 ) 


1.137 
0.971 
0.955 


Standard fit 


1607 1772 2182 2841 


( 1365\ 


(\ 7 7°4 8 ) 


fl73l\ 

V 136 / 


(2in\ 
l 345 J 


0.958 


Sol a = 0.1 
Sol a = 0.5 
Sol a = 1.0 


1606 1772 2182 2796 

1607 1772 2184 2811 
1607 1772 2186 2784 


( 1367\ 

( 1369\ 
V 173 j 
/ 1374\ 


(Z) 


( 1732\ 
V 13 6y ) 
( 1730\ 

V 135,1 
( 1730\ 

V "5 ) 


/2115\ 
^ 338 ) 
f 212l\ 
1, 343 ,1 
/2127\ 
^ 362 J 


1.001 
0.973 
0.995 



range 2 GeV 2 < s < 5 GeV 2 . 
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TABLE VII: The extracted D13 partial wave T-matrix poles using different solutions of channel propagator. We change interme- 
diate part of channel propagator and the modification is the same for all three used channels 



Solutions 


Rorp TiolpQ 
±J Cl.l C IJU1CB 


Dressed poles 


Xr 


W S1 W S2 W S3 
MeV 


(■ 


RcW \ f 
-2ImW J \- 


RcW \ ( 
-2ImW ) \- 

MeV 


RcW \ ( 
-2ImwJ y." 


ReW \ 
-2ImW ) 


Sol a = -1.0 
Sol a = -0.5 
Sol a = -0.1 


1578 1881 2503 

1579 1881 2500 
1581 1880 2498 




fl502> 

( 150^ 
V 113; 
fl505> 
V 120y 




( 1806> 
V123; 
f 1807N 
V 125y 
fl808^ 




fl937^ 
v 384y 

y 436 y 
^194lA 

V 479 y 




( 2675^ 
^ 583 y 

(imi\ 

v 567 y 
f 2697^ 
v 575 y 




1.396 
1.087 
1.027 


Standard fit 


1582 1880 2499 




fl506^ 
V 121y 




( 1807A 
V127y 


(IT) 


^269^ 
y 583 y 




1.027 


Sol a = 0.1 
Sol a = 0.5 
Sol a = 1.0 


1582 1879 2497 

1584 1878 2495 

1585 1876 2493 




^1505^ 
V 121y 
fl506^ 
V127y 
( 1507N 
^135; 




^1808^ 
V127y 
fl808^ 
^127y 
(\%\\\ 
V127y 




f 1939N 
v 500 y 
^1935^1 
v 546 y 
? 1929A 
V 610y 




^2685^ 
v 594 y 
f 2670A 
v 609 y 
f 2656N 
V 627y 




1.032 
1.084 
1.183 



500 
400 
' 300 
200 
100 



1 1 1 1 

: (a) 


1 1 1 1 


- 


"■1,1, 
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1850 



400 
350 
I 300 
^250 
J 200 
150 
I0i 
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FIG. 7: (Color online) The extracted Sn , P\\ and D13 partial waves T-matrix poles using different solutions of channel propagator. 
We change intermediate part of channel propagator and the modification is the same for all three used channels. Red arrows 
indicate the position of Standard solution poles. 
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We observe: 

1. The strong modification of the imaginary part of the channel propagator is via the dispersion relation 



Eq. ( 10 ) transmitted to the real part, causing opposite effect (reduction) for the S-wave, and a combined 



effect for other partial waves. 

2. From Tables |V} | VII| we learn that bare poles are fairly stable with respect with change of the intermediate 
part of the channel propagator, while dressed poles experience quite a change. 

3. The dressed poles can, however, be classified in two categories: 

a. Poles weekly dependent upon the strong changes of the intermediate part - stable poles 

Sn (1518), Su (1642), Pn (1731), D 13 (1506) and D 13 (1807) 

b. Poles moderately and notably dependent upon the strong changes of the intermediate part - running 

poles 

Su (1779), Pu(1365), P u (1731), P u (2117), Di 3 (1939) and D 13 (2691) 

4. We may connect the category of the pole to its self energy content: 

a. For the stable poles the self energy correction to the bare mass is small 

b. For the running poles, the self energy correction to the bare mass is either strong (as for Su (1779), 

Pu (1731), Pu (2117) and Di 3 (2691) ), or poles are entirely dynamically generated meaning that 
they do not at all have a corresponding, nearby bare poles (as for Pu(1365) Roper resonance, and 
Di 3 (1939) ) 

5. Consequently, "difference" of bare and dressed pole position can be taken as the criterion for the stability 
of the pole position with respect to the strong modification of the intermediate part: "The lower the self 
energy correction is, the more stable the pole position is." 

6. As seen in Tables |v| - |VIl| the reduced Xr> the quantity representing the measure of ability of our model 
to reproduce the input data, is systematically smallest for the non-modified solution. In other words, 
our fit to the input data worsens whenever any change to the intermediate part of channel propagator is 
introduced. 



3. Threshold behavior 



In this part of our analysis, behavior of channel propagator's imaginary part close to the threshold a given 
channel is changed. For values s > 3.5 GeV 2 channel propagator is kept unchanged. The threshold behavior of 
imaginary part of the channel propagator is given by q 2L+1 . It is changed by multiplying the channel propagator 
with q x . The continuity in energy is ensured by the matching function h(s) = (a + b ■ s + c ■ s 2 ). Modified 
imaginary part is given by: 

Im$corr = Im$ • • h(s). (20) 

with A = ±0.5. Parameters a, b and c are adjusted in such a way to ensure continuity of channel propagator 
and first derivative at s =3.5 GeV 2 . 

The change of the channel propagator for the Su, Pu and Di 3 partial wave is given in Fig. ([8]). 
Result: 



The shift in pole positions for partial waves Su, Pu and Di 3 are presented in Tables VIII Xjand depicted 
in Fig. ([9]). We personally do not believe that such strong variations in the threshold behavior are physics- 
wise permissible, but we have investigated this aspect of our model just to check the how quickly our results 
deteriorate when even uphysical changes are introduced. We are happy to conclude that the results are only 
moderately sensitive to even huge changes. 
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FIG. 8: (Color online) Real and imaginary parts of the Sn, Pn and D13 elastic channel propagator. Threshold behavior. Sol 
thl: Im<E> corr = Im<!> ■ yfq{a + cs 2 ), Sol th2: Im<l> corr = Im<3? ■ ^/q{a + 0.5s + cs 2 ), Sol th3: Im3> corr = Im<I> ■ (-^=)(a + cs 2 ), Sol th4: 

Im<l> corr = Im«l> • (-^=)(a — 0.5s + cs 2 ). 

TABLE VIII: The extracted Sn partial wave T-matrix poles using different solutions of channel propagator. We change threshold 
behavior of channel propagator and the modification is the same for all three used channels. 



Solutions 



Bare poles 



W Sl W S2 W s 
MeV 



Dressed poles 



/ RcW \ 


( RcW \ 


( RcW \ 


I -2ImW J 


I -2ImW J 


I -2ImW J 



MeV 



Xr 



Standard fit 



1523 1640 1834 



1518 
189 



1642 
208 



1779 
427 



1.016 



Sol thl 
Sol th2 
Sol th3 
Sol th4 



1551 1673 1899 
1557 1683 1915 
1477 1575 1728 
1453 1540 1678 



1505 

139 
1505 

130 
1552 : 

164 
1557 

154 



1642 
212 
1642 
208 
1642 
212 
1643 
203 



1755 
468 
1751 
457 
1781 
397 
1786 
401 



1.289 
1.372 
1.409 
1.587 



TABLE IX: The extracted Pn partial wave T-matrix poles using different solutions of channel propagator. We change threshold 
behavior of channel propagator and the modification is the same for all three used channels 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W B3 W S4 
MeV 


/ ReW \ 
^-2ImW/ 


f ReW \ f 
\-2JmWj V 
MeV 


RcW \ 
-2ImW J 


( ReW \ 

\-ii-am) 


Xr 


Standard fit 


1607 1772 2182 2841 


/ 1365\ 
{l57j 


Cm) 


>T73l\ 

v i 36 ; 


(2in\ 

V 345 J 


0.957 


Sol thl 
Sol th2 
Sol th3 
Sol th4 


1629 1780 2192 2802 
1639 1784 2196 2801 
1572 1763 2174 2712 
1534 1753 2168 2546 


/T376\ 
^ 161 ) 
( 1383\ 

v ^3 ; 

/l35l\ 

V 163 ) 

('.«) 


(Z) 

(l64 6 ) 
(-) 

C 7 9 \ 8 ) 


^1733\ 
y 12SJ 
^1734\ 
V 12 4y ) 
^1729\ 
V 14 3 J 
^1724\ 
V 145 / 


/2112\ 
^ 327 ) 
/2112\ 
^ 323 ) 
/2123\ 
^ 352 ) 

/2135\ 
V 344J 


1.157 
1.335 
1.010 
1.274 



TABLE X: The extracted D13 partial wave T-matrix poles using different solution of channel propagator. We change threshold 
behavior of channel propagator and the modification is the same for all three used channels 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W S3 
MeV 


/ ReW \ 
l^-2ImW J 


( ReW \ f 
^-2ImW J \- 

MeV 


RoW \ 
-2ImW J 


f ReW \ 
l v -2ImWy 


Xr 


Standard fit 


1582 1880 2499 


A506\ 
{l2l) 




' 1939\ 
v 48 5 y ) 


/269l\ 
^ 583 J 


1.027 


Sol thl 
Sol th2 
Sol th3 
Sol th.4 


1589 1881 2500 
1593 1882 2505 
1573 1875 2481 
1563 1871 2468 


A503\ 
1,114 ) 
/1502\ 
{l09j 
/1508\ 

\ 12S ) 
A515\ 
{137) 


CEO 
Ci 8 2 ° 9 8 ) 


'1933\ 

V 449 / 
f 193u\ 
, 424,1 

f 1947\ 
,518,1 
/ 196l\ 
v 55 4 y ) 


/2694\ 
{ 567 ) 
/2692\ 
1, 563 ; 
/2680\ 
1^609 J 
f2662\ 
1^648 J 


1.057 
1.148 
1.119 
1.363 
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FIG. 9: (Color online) The extracted Sn, Pn and D13 partial waves T-matrix poles using different solutions of channel propagator. 
We change threshold behavior of channel propagator and the modification is the same for all three used channels. Red arrows 
indicate the position of Standard solution poles. 
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Conclusion: 



As seen in Fig. (|8) the threshold behavior was modified by the continuous and arbitrary function given by 
Eq. Q. We have c 



losen 4 possibilities, in two of them the threshold behavior has been made shallower (thl 
and th2), and for the two of them (th3 and th4) much steeper. 

We observe: 

1. The strong modification of the imaginary part threshold behavior of the channel propagator is via the 



dispersion relation Eq. (10) transmitted to the real part, causing significant changes of the real part in 
complete energy range. 

2. Contrary to the former case when intermediate part of the channel propagator has been analyzed, bare 
poles now experience significant change for all energies, and all partial waves. 

3. Dressing tries to compensate for the strong shift of the bare poles, but not always very successfully. 

4. As a result, we conclude that the correct threshold behavior of the channel propagator is essential for the 
applicability of the model. 

5. With some possible coincidental exceptions, all dressed poles experience a notable shift. 



6. As seen in Tables VIII - DO the reduced Xr, the quantity representing the measure of ability of our model 
to reproduce the input data, is systematically smallest for the non-modified solution. In other words, 
our fit to the input data worsen whenever any change to the threshold behavior of channel propagator is 
introduced. 



C. Dispersion integral subtraction constant 



As positions of bare poles are in principle not too well defined (they must depend on the renormalization 
scheme, see Capstick et.al the fact that they can only be significantly moved by changing the threshold 

behavior troubled us a lot. As our model is an effective model, changing its ingredients should enable us to 
"mimic" at least some of other models, and that means we should observe a significant shift of bare poles. So, 
having the bare pole model dependence hidden only in fairly well defined threshold behavior was non-satisfactory. 

The answer to this issue was already hinted by Vrana et al in [T3] where they have shown that position 
of bare poles significantly depends on the dispersion integral subtraction constant, and its importance for 
the position of the bare and dressed poles was briefly discussed. We have, therefore, extensively tested this 



freedom within our model, and shifted the subtraction constant Re$(s ) defined in Eq. (10) from its original 
zero value. We have chosen SC Sol 1: Re$(s ) = for all three used channels; SC Sol 2: Re<I>(so) = for 
irN channel only, while in other two channels Re$(so) = 0; SC Sol 3: Re<I>(so) = m n f° r channel only, 
while in other two channels Re$(so) = 0; SC Sol 4: Re<I>(so) = for ir 2 N channel only, while in other two 
channels Re$(so) = 0; SC Sol 5: Re$(so) = — for all three used channels, SC Sol 6: Re<&(so) = fn a , where 
m a is the mass of the meson in corresponding channel. 

As changing the dispersion integral subtraction constant only shifts the real part of the channel propagator 
into positive or negative y-direction for a certain constant value, we show the modified real part of the Sn, Pn 



and D13 partial wave channel propagator in Fig. ( 10 1 for only two representative values of subtraction constant 
Re<I>(so) = irn-Tr in all three channels simultaneously. For other modifications related plots look qualitatively 
the same, so we do not show them explicitly. 

Results: 



Results are given in Tables XI -XIII and in Fig. ( 11 1. 



We see that bare pole positions are notably changed with the variation of the dispersion relation sub- 
traction point value. For the shift into the time-like direction the bare poles are shifted upwards; and 
when the subtraction point is shifted into the negative energy region the value of bare poles is reduced. 
The change is notable, and for inspected values varies for 150 - 200 MeV. It is important to notice that the 



dressed pole positions are extremely stable, see Fig. ( 11 ), and not influenced by the subtraction constant values. 
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FIG. 10: (Color online) Real part of Six, P-i an( l f 13 elastic channel propagator. We change the dispersion integral subtraction 
constant. 
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FIG. 11: (Color online) The extracted Sn , Pn andDi3 partial waves T-matrix poles using different solutions of channel propagator. 
We change dispersion integral subtraction constant for all three used channels. Red arrows indicate the position of Standard solution 
poles. 
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TABLE XI: The extracted Sn partial wave T-matrix poles using different values of subtraction constant. SC Sol 1: Re$(so) = m-ir 
for all three used channels; SC Sol 2: Re<J>(so) = for nN channel only, while in other two channels Re<I>(so) = 0; SC Sol 3: 
Rc<J>(so) = m v for r)N channel only, while in other two channels Re<I>(so) = 0; SC Sol 4: Rc<I > (so) = m ^ f° r n 2 N channel only, 
while in other two channels Rc<I>(so) = 0; SC Sol 5: Re<I>(so) = — for all three used channels, SC Sol 6: Re$(so) = m a , where 
m a is the mass of the meson in corresponding channel. 





Bare poles 


Dressed poles 






W S1 


w 32 

MeV 


w S3 


( RcW 
l,-2ImW 


\ ( RcW \ ( 
j l,-2ImW J 1 

MeV 


' ReW \ 
v -2ImW; 


Xr 


Standard fit 


1523 


1640 


1834 


( 1518\ 

V 189 / 


f 1642\ 
1, 208 ; 


/l779\ 
1, 427 ; 


1.016 


SC Sol 1 
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1686 


1953 


f 1520\ 
I 188 J 


/ 1643\ 

V 216 ; 


/l776\ 
1,456 ; 


1.054 


SC Sol 2 


1537 


1662 


1868 


/1520\ 

1 189; 


/1642\ 
1,210 J 


/1778\ 

I 433 ; 


1.019 


SC Sol 3 


1587 


1730 


1907 


( 1518\ 
V188y< 


/ 1642\ 
1, 207 ; 


/l779\ 
1, 425 ; 


1.016 


SC Sol 4 


1529 


1659 


2150 


( 1522\ 

ViseJ 


/ 1647\ 
{216) 


/ 1774\ 

v 493 ; 


1.141 


SC Sol 5 


1462 


1593 


1756 


(Sf) 


( 1643\ 


/1776\ 

^406 ; 


1.052 


SC Sol 6 


1663 


1761 


2252 


( 1523\ 
V 186^ 


( 1649\ 

V2i8; 


/ 1775 \ 

v 510 ; 


1.216 



TABLE XII: The extracted P\\ partial wave T-matrix poles using different values of subtraction constant. SC Sol 1: Re<I>(so) = win 
for all three used channels; SC Sol 2: Re<J>(so) = win for nN channel only, while in other two channels Re$(so) = 0; SC Sol 3: 
Re<J>(so) = ™>7 for r)N channel only, while in other two channels Rc<3>(so) = 0; SC Sol 4: Rc<I>(so) = m ^ f° r n 2 N channel only, 
while in other two channels Re<I>(so) = 0; SC Sol 5: Re«I>(so) = — for all three used channels, SC Sol 6: Re<J>(so) = m a , where 
m a is the mass of the meson in corresponding channel. 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W S3 W S4 
MeV 


/ ReW \ 
l v -2ImW; 


/ ReW \ / 
^-2ImW; 1,- 

MeV 


RcW \ 
-2ImW ; 


/ ReW \ 
^-2ImW; 


Xr 
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0.965 


SC Sol 3 


1611 1796 2376 3586 


/ 1364\ 
V159; 


(-) 


{ 1730\ 
V13 5 ; 


/2116\ 
1, 349 ; 
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0.959 


SC Sol 5 


1461 1634 1966 2412 


/ 1365\ 
\W5) 


(Ye") 


A73l\ 

V i 40 ; 


/2110\ 

1, 325 ; 
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1, 360 ; 
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TABLE XIII: The extracted D13 partial wave T-matrix poles using different values of subtraction constant. SC Sol 1: Re<I>(so) = 
m-n for all three used channels; SC Sol 2: Re<S>(so) = "V for channel only, while in other two channels Rc<I>(so) = 0; SC Sol 3: 
Re<J>(so) = m v for r]N channel only, while in other two channels Rc<I>(so) = 0; SC Sol 4: Re<I>(so) = "V2 for tt 2 N channel only, 
while in other two channels Re<I>(so) = 0; SC Sol 5: Re<I>(so) = — m n for all three used channels, SC Sol 6: Re$(so) = m a , where 
m a is the mass of the meson in corresponding channel. 





Bare poles 


Dressed poles 
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We have also tested the freedom in choosing the dispersion integral subtraction point point sq. In Fig. (12) 
where the resulting real part is depicted, we show that the results are qualitatively and quantitatively similar to 
changing the size of the subtraction constant. We see that changing the subtraction point value shifts the real 
part upwards and downwards depending on the direction of the shift. So, in practice, this freedom is equivalent 
to retaining the subtraction point at threshold and varying the subtraction constant. Therefore, the similar 
conclusions as in case of dispersion constant variation are drawn that the bare poles are strongly influenced, 
and dressed poles remain virtually untouched. 
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FIG. 12: (Color online) Real part of Sn, Pxi and D13 elastic channel propagator. We change the dispersion integral subtraction 
point value. 



Conclusion: 

Results shown in Tables |Xl] - pOTT] prove that the value of subtraction constant affect bare pole positions 
significantly. Dressed poles remain practically unchanged. 



D. Background parameterization 



Additional background poles below and above physical region are added. Our model was based on the premises 
that two subthreshold unphysical poles are enough to represent the smooth energy dependent background, and 
we test this assumption. We test two aspects: 

1. Are two poles absolutely essential? 

2. Does adding new subthreshold, or distant, positive energy poles change the pole positions? 

Therefore, we have added first one, and then two poles, below, and far above threshold, and with different 
signs representing attraction or repulsion. All possible combinations for three analyzed partial waves are 
analyzed. 

Result: 



The shift in pole positions for partial waves Sn, Pn and D 13 are presented in Tables |XIV] - |XVI| and depicted 
in Fig. (pi. 
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TABLE XIV: The extracted Sn partial wave T-matrix poles. Additional background poles below and above physical region are 
added. 







Bare poles 




Dressed poles 






Solutions 


Bl 
B2 
GeV 2 


B3 
B4 

GeV 2 


W B1 W S2 W S3 
MeV 


/ ReW 
I -2ImW 


\ t RcW \ / 
j \ -2ImW / I 

MeV 


RoW \ 
-2ImW J 


Xr 


Standard fit 


-1.256 
0.210 


1523 1640 1834 


/1518\ 

V 189 ) 


/1642\ 
I 208 J 


/1779\ 
V 427 7 


1.016 


el = +\(Below) 

jui J. e^ — — j. i D niu w j 
e3 = +l{Below) 


-1.382 
0.460 


-1.851 
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\™») 
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/1779\ 
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l ni is 
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Snl 9 p9 — —1 ? RrlmiA 

kJUl Zi — A.yUCHJUJ) 

e3 — — l(Se/oiu) 


-2.508 
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1523 1640 1834 
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SQ1 J e3 = -l(Beioiu) 
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e2 — — l(Above) 
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el = +\{Below) 
Sol 6 e2 = -l(Above) 
e3 = +l(Afco«e) 


-7725 
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TABLE XV: The extracted P\\ partial wave T-matrix poles. Additional background poles below and above physical region are 
added. 



Solutions 



Bare poles 



Bl B3 
B2 B4 
GeV 2 GeV 2 



W S1 W S2 W S3 W S4 
MeV 



Dressed poles 



( RcW \ 


/ ReW \ 


( RoW \ 


/ ReW \ 


I -2ImW) 


I -2ImW j 


I -2ImW) 


I -2ImW J 



MeV 
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TABLE XVI: The extracted D13 partial wave T-matrix poles. Additional background poles below and above physical region are 
added. 



Solutions 



Bare poles 



Bl 
B2 

GeV 2 



B3 
B4 
GeV 2 



W S1 W S2 Ws 3 
MeV 
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-2ImW 



RcW \ ( RcW 
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MeV 
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FIG. 13: (Color online) The extracted 5n, P\i and D13 partial waves T-matrix poles. Additional background poles below and 
above physical region are added. Red arrows indicate the position of Standard solution poles. 
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Conclusion: 



XIV 


XVI 


in Fig. ( 
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(131 we see that all extracted bare and dressed poles are practically iden- 
tical with the exception of solutions Sol 4 and Sol 6. And these are solutions where only 1 subthreshold pole 
has been allowed. 



Therefore: 

1. at least two background poles are needed to adequately represent the energy dependent background. 

2. two subthreshold poles are enough to adequately represent the background. Two pole meromorphic 
approximation for the background representation is justified. 



E. Number of channels 



It is not self evident that introducing new channels, i.e. new contributions to the self energy term in the 
denominator of Eq. ( 11 1, does not influence pole positions at all. Namely, self energy matrix is of higher order, 



and self energy is distributed in a different way. To test this assumption we have introduced the fourth channel, 
and fitted the existing data base with this channel. 



In addition to the existing irN elastic and irN — » r/N processes we have introduced a third channel opening 
at 1612 MeV (KA channel). The mass of the effective channel has not been altered, as the effective channel is 
primarily dominated by the 7T7tN channel. The results are given in Tables XVII - XIX 



TABLE XVII: Extracted Sn partial wave T- matrix poles obtained by including three and four channels into a model. 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W S3 
MeV 


( RcW 
l v -2ImW 


\ f RcW \ f 
) \-2lmW) \ 

MeV 


RoW \ 
-2ImW ) 


xl 


Sol 3CH 


1523 1640 1834 


/1518\ 
V 189 J 


/l642\ 
^ 208 ) 


( 1779\ 
{427) 


1.016 


Sol 4CH 


1523 1640 1836 


/1516\ 
{l90j 


/1642\ 


/1776\ 
{418 ) 


1.013 



TABLE XVIII: Extracted P\\ partial wave T-matrix poles obtained by including three and four channels into a model. 





Bare poles 


Dressed poles 




Solutions 


W S1 W S2 W S3 W S4 
MeV 


( RcW \ 
\-2lnw) 


/ ReW \ f RcW \ 
l v -2ImWj l v -2ImWy 

MeV 


( RcW \ 
l v -2ImW/ 


Xr 


Sol 3CH 


1607 1772 2182 2841 


/l365\ 
1 157 ) 


/1708\ /1731\ 
^ 174 J \ 136 J 


/2117\ 
I 345 ) 


0.958 


Sol 4CH 


1607 1772 2182 2841 


/1363\ 
1,161 J 


/1708\ /1731\ 
^ 174 J \ 136 J 


(2in\ 

V343j 


0.970 



Conclusion: 

Pole position for the analyzed partial waves are extremely stable with respect to the opening of a new 
channel. 



TABLE XIX: Extracted D13 partial wave T-matrix poles obtained by including three and four channels into a model. 





Bare poles 


Dressed poles 




Solutions 


w 81 W S2 W S3 
MeV 


( RcW \ 
\-2\mW) 


( RcW \ ( RoW \ 
l v -2ImW ) l v -2ImW J 

MeV 


( RoW \ 
l v -2ImW J 


Xr 


Sol 3CH 


1582 1880 2499 


( 1506\ 
\l2l) 


/ 1807\ / 1939\ 
^ 127 J \ 485 J 


/269l\ 
^ 583 J 


1.027 


Sol 4CH 


1581 1890 2499 


/1500\ 
\ 169 ) 


A808\ /1944\ 
^ 131 J \ 492 J 


/2689\ 
^ 586 J 


1.035 
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F. Mass of the effective channel 



However, in spite of looking identical to all other tests, choice of the effective cannel mass undergoes strong 
physical constraints: 

1. The effective channel represents all inelastic two and three body channels as only one, stable, two body 
channel. This is a serious assumption because some channels, the dominant three body channel like ttt:N 
and unstable ones like aN and pN introduce a significantly different analytic structure, so this might 
influence the conclusion on existence of some resonant states. An indication that it might be actually 
happening is already given for the pN channel where its complex branchpoint for the Pn partial wave 
lies dangerously close to the Pn(1710) resonance - see ref. [53]. Genuine three body effects for the 
dominant irirN channel might also cause serious problem, This particularities should be analyzed as they 
are detected, but the method in general can enable statistical analyses of most of other resonant states. 

2. As the effective channel represents the loss of flux to all inelastic channels, its mass has to be neither to 
high nor too low. If it is too high, the loss of flux to inelastic channels imposed by the data at energies 
below effective channel threshold will not be compensated by the third channel, so the T-matrix values 
will simply be too high to fit the data. If the mass is too low, we shall artificially open the new degrees 
of freedom where they actually do not exist. So, the mass of the effective channel will be fairly precisely 
defined by all inelastic channels. 

3. When the effective channel mass is changed, the effective channel T-matrix threshold is also shifted. 
Therefore, the code will unsuccessfully try to reproduce the inelastic channel data which are generated by 
the different threshold, and move other resonances. 

Therefore, we have decided to show both aspects of the problem for the Sn, P\\ and D13 partial waves: In 



Fig ( 14 ) and Tables XX - XXII we illustrate what is happening to the poles when the mass of the effective 
channel is changed ± 40%. 




FIG. 14: (Color online) Real and imaginary parts of S\\, P11 and D13 effective channel propagator. 



Result: 



The shift in pole positions for partial waves Sn, Pn and D13 are presented in Tables |XXf |XXII| and 
depicted in Fig. (15 1. 
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TABLE XX: The extracted Sn partial wave T-matrix poles using different mass of effective channel 



Solutions 



Bare poles 



W S1 W S2 W B3 
MeV 



Dressed poles 



RcW \ ( RcW \ 
-llmW) \-2lmW) 

MeV 



ReW 
2ImW 



Xr 



Sol 0.6mEF 
Sol 0.7m EF 
Sol 0.8m BF 
Sol 0.9m EF 



1522 1637 1823 

1523 1638 1825 
1523 1638 1827 
1523 1639 1830 



1520 
184 

1520 
184 

1520 
185 

1519 
186 



1644 
210 
1644 
209 
1644 
210 
1643 
209 



1776 
430 
1777 
430 
1776 
429 
1779 
426 



1.028 
1.024 
1.020 
1.014 



Standard fit 



1523 1640 1834 



1518 
189 



1642 
208 



1779 
427 



1.016 



Sol l.lm_BF 
Sol 1.2m_BF 
Sol 1.3m,EF 
Sol lArriEF 



1524 1641 1837 
1524 1643 1842 
1523 1643 1846 
1522 1644 1850 



1518 
188 
1516 
188 
1511 
189 
/ 1499 
^ 185 



1642 
207 
1642 
206 
1640 
202 
1638 
199 



\ / 1638\ / 179l\ 
J \ 199 J \ 445 J 



1778 
427 
1781 
427 
1786 
437 
1791 
445 



1.003 
1.016 
1.062 
1.172 



TABLE XXI: The extracted Pn partial wave T-matrix poles using different mass of the effective channel. 



Solutions 



Bare poles 



W S1 W S2 W S3 W S4 
MeV 



Dressed poles 



RcW 
-2ImW 



RcW \ ( RcW 
-2ImWj l v -2ImW 

MeV 



ReW \ 
-2ImW I 



XR 



Sol 0.6m,EF 
Sol 0.7m,EF 
Sol 0.8m EF 
Sol 0.9m E F 



1574 1775 
1582 1774 
1595 1772 
1600 1772 



2204 2826 
2199 2811 
2182 2773 
2178 2799 



1348 
220 
1347 
215 
1344 
216 
1353 
192 



1725 

235 
1722 

218 
1718 

183 
1714 

169 



1729 

148 
1728 

144 
1725 

145 
1729 

145 



2131 

383 
2129 

376 
2130 

347 
2115 

344 



1.827 
1.587 
1.272 
1.002 



Standard fit 



1607 1772 2182 2841 



1365 
157 



1708 
174 



1731 
136 



2117 
345 



0.958 



Sol 1.1m ef 
Sol \.2m E F 
Sol 1.3m_BF 
Sol lArriEF 



1614 1773 

1615 1775 
1608 1777 
1592 1779 



2184 2811 
2182 2809 
2175 2832 
2166 2883 



1377 

132 
1399 

120 
1423 

114 
1447 

119 



1707 

163 
1705 

156 
1704 

158 
1702 

162 



/1447\ /1702\ /1747\ 

^ 119 J \ 162 ) \ 137 ) 



1732 
132 

1735 
132 

1735 
132 
1747 
137 



2109 
328 
2103\ 
320 ) 
2103 
320 
( 2092\ 
^ 330 ) 



1.182 
2.215 
4.166 
7.021 



Conclusion: 

We see that increasing the mass of the effective channel too much above the threshold of the first inelastic 
channel contributing to the particular partial wave is producing unrealistic shifts of the pole positions. 

Reducing the mass for up to 40% is producing some shift, but not significant. And again we see that the 
T-matrix poles form two classes: the first class-almost stable poles Sn (1518), Sn (1642), D13 (1506), D13 
(1807), and D 13 (2691) and the second class- more sensitive poles Sn (1779), Pu (1506), Pn (1708), P n 
(1731), Pn (2117), and D 13 (1939). The pole is not sensitive to the model parameters if it is dominantly 
determined by the bare singularity (small self energy correction), and it is moderately or strongly sensitive 
upon the model choices if the pole is created more dynamically (self energy contribution is strong) . 
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TABLE XXII: The extracted D13 partial wave T-matrix poles using different mass of the effective channel 



Solutions 


Bare poles 


Dressed poles 


xl 


W S1 W S2 W S3 

MpV 


(- 


ReW \ / 
-2lmW) \- 


ReW \ / 
-2lmWJ \- 

MeV 


ReW \ f 
-2ImWy l l v - 


RcW \ 
-2ImWy l 


Sol O.QjTlEF 

Sol 0.7m E F 
Sol 0.8m EF 
Sol 0.9m BF 


1580 1878 2496 
1580 1878 2496 
1580 1879 2499 
1582 1880 2499 




( 1509^ 
V129y 
(\bm\ 
I 129 ; 

( 1507> 
V 127y 

(\bm\ 




V 129y 
^1808^ 
I 128 ; 

V128y 
(1807\ 
V 127y 




f 1945^ 
^ 526 y 

\ 520 / 

^194^ 

V 512y 
^ 194CA 
^ 488 y 




^2687^ 
^ 588 y 
f2687> 

V 589 y 
f2688> 
^ 589 y 

/^gi^i 

V 578 y 




1.096 
1.069 
1.041 
1.044 


Standard fit 


1582 1880 2499 




( 1506> 
V 121y 




^1807^ 
V 127y 




f 1939^ 
V 485y 




^269^ 
V_ 583 y 




1.027 


Sol 1.1m ef 
Sol 1.2m_BF 
Sol 1.3iriEF 
Sol IAulef 


1583 1880 2499 

1584 1881 2501 
1586 1882 2504 
1589 1883 2507 




f 1505^ 
V118y 

V 113y 
(lb02\ 

(lbOl\ 




V 127y 
^1806^ 
V 126y 

V 123y 
( 1804> 

V 12 °y 




( 1936^ 

V 472 ) 

( 1935A 

V 452 y 
f 1935A 

V 426y 
( 1934A 

V 397 ) 




^2694^ 

V 572 y 
^2690A 

V 570 y 

V 571y 
f 267CA 
^ 573y 




0.996 
1.031 
1.127 
1.325 



500 
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100 



1 1 1 1 1 1 1 1 1 1 


1 : 


: i i 
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FIG. 15: (Color online) The extracted Sn, P\\ and D13 partial waves T-matrix poles using different mass of the effective channel. 
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From Tables XX|XXII it is seen that reduced \ 2 significantly raises when we increase the mass of the effective 
channel for up to 40%. This is, however, not so when the mass of effective channel is decreased for the same 
amount. We can easily understand the effect. Namely, the effective channel serves for ensuring unitarity for 
all open inelastic channels. However, if we increase the mass of the inelastic channel too much, then the fit, 
which is manifestly unitary, will not have an extra channel at its disposal to get rid of unwanted flux, and the 
code will be forced to fit the input data with amplitudes which are much more "elastic" then in reality. Simply, 
there will be no room for the loss of flux to inelastic channels which are open at these energies. Therefore, the 
fit will produce unrealistic results. On the other hand, if we lower the mass of the effective channel for 40%, 
we shall allow the fit to send the flux to inelastic channel if needed. However, the fit declines this possibility 
because the input data are dominantly elastic. 



VI. CONCLUSIONS 



In our model we find two classes of poles: ones strongly determined by the bare poles, and ones which are 
more influenced by the self energy correction. 

Poles Sn (1518), Sn (1642), P n (1731), D 13 (1506) and D i3 (1807) belong to the first class, and Sn (1779), 
Pn(1365), Pn (1708), Pn (2117), D 13 (1939) and D 13 (2691) to the second. 

We show that the position of bare poles is model dependent. It is strongly influenced by dispersion relation 
subtraction constant and the channel propagator threshold behavior. As we believe that the threshold behavior 
is fairly well determined by the q 2L+1 law, the major model freedom of bare poles is originating from the choice 
of dispersion relation subtraction constant. 

The dressed pole positions are much better fixed. Poles for which the self energy correction is small are 
stable, and ones for which the self energy correction is stronger are only moderately running. The "complex 
plane distance" between bare and dressed poles representing the size of the self energy contributions can be 
taken as a measure how model dependent a certain pole is. The smaller the difference is, the more confident 
the pole extraction is. 



We observe that the reduced Xr-> the quantity representing the measure of ability of our model to reproduce 
the input data, is systematically smallest for the non-modified solution. In other words, our model is unable 
to correctly fit the input data whenever any significant change to the model ingredients is introduced. We 
conclude that our model is sensitive to the analytic structure embedded into the input data by the way they 
are generated, because the underlying analytical structure of the Zagreb model is transferred to the data when 
we created the input data set by distributing Zagreb theoretical curves. If our analytic structure is wrong, and 
it certainly has to be up to a certain level, the real data coming from the experiment are described by the 
different analytic function and will not be ideally fitted by our curves. So, pole positions will be shifted. 

However, we could at least in principle, use this feature as to test the shape of the model ingredients deter- 
mining the self energy corrections. We propose the following procedure to improve the analytic structure of 
Zagreb CMB model with respect to any input data, either originating from different PWA or preferably coming 
directly from experiment in the form of partial wave data: 

1. We fit the obtained input set and get the pole positions with the standard channel propagator, and note 
the belonging \r- 

2. We change the form of the channel propagator, re do the fit, and compare the new \r with the old value. 

3. If the Xr is better, we keep on changing the form of the channel propagator because we infer that our 
imposed analytic structure is not adequate, and that we are finding the better one. 

4. When we achieve the best Xr> we pronounce that we have found the optimal value of the channel propa- 
gator for our formalism, and declare the found pole positions a correct ones. 

5. We should, in principle, be able to follow which poles are changed, and which remain stable, and on the 
basis of that infer the error. 

However, the afore hinted procedure is not within the scope of this paper. 
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In conclusion: 

Zagreb CMB model can be used for extracting dressed pole positions from either partial wave data 
or partial wave amplitudes. The bare pole masses are strongly model dependent, while the dressed pole 
positions are much better fixed. The poles which are strongly influenced by the self energy corrections, 
however, show more model dependence then the ones which are dominated by the bare poles. At this point 
it is important to mention that our findings about the relative size of bare vs. self energy contributions of a 
certain pole are not eternal truth for all models. As bare mass is a movable quantity, it is easy to imagine 
that for some other model the shift of the bare pole position in proper direction might result in reducing the 
importance of self energy corrections. So, poles which are dominated by self energy corrections in our model 
might be dominated by the bare pole in others. This allows us to speculate that changing the subtraction 
constant that is a completely free parameter of our model may shift the results of Zagreb CMB fit in the 
direction of other models with respect to the value of bare poles. 

There are other interesting questions like the role of physical requirements on the S-matrix that are missing 
in the Zagreb CMB model. Those are, e.g., sub-threshold cuts (circular, short nucleon, left-hand cuts as we 
have a trivial left-hand cut from kinematics, but that is, of course, not the full one...), or true and not effective 
three-body intermediate states with the corresponding additional analytic structures in the complex plane. 
Those are the really questions of exceptional interest which are really hard to answer, and will not be answered 
yet in this paper. 
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